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In the past few decades, several models have predicted an energy-dependence of the speed 
of light in the context of quantum gravity. For cosmological sources such as active galaxies, 
this minuscule effect can add up to measurable photon-energy dependent time lags. In this 
paper a search for such time lags during the H.E.S.S. observations of the exceptional very 
high energy flare of the active galaxy PKS 2155-304 on 28 July in 2006 is presented. Since no 
significant time lag is found, lower limits on the energy scale of speed of light modifications 
are derived. 

PACS numbers: 12.60.Jv, 04.60.-m,11.25.Wx,96.50.S- 



Albert Einstein's postulate "that light is always prop- 
agated in empty space with a definite velocity c which is 
independent of the state of motion of the emitting body" 
[H is one of the pillars of modern physics. Modification 
of this postulate would have far-reaching consequences 
for our understanding of nature, it is therefore impor- 
tant to constantly improve the verification of its validity. 
Particularly in the past few decades, a possible energy 
dependence of the speed of light has been predicted in 
the framework of quantum gravity models [3, [1, 01 and 
effective field theory Q, leading to deviations from this 
postulate (for reviews see [^, The speed of light 

modifications have different functional dependencies on 
the photon energy and helicity in different models. Pre- 
dictions usually entail free parameters such as the rele- 
vant mass scale. However, it is commonly expected that 
this modification should appear at energies of the order of 
the Planck energy {Ep = 1.22 x 10^^ GeV). For energies 
much smaller than the Planck energy, a series expansion 
is therefore expected to be applicable, allowing the en- 
ergy dependence of the speed of light to be parameterized 
in a model-independent wayQ. The photon speed c' is 
written up to second order in energy E as: 



1 + + c— 
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where ^ and C are free parameters. Even for the highest 
photon energies currently measured the corrections are 
expected to be very small. However, Amelino-Camelia 
et al. @ suggested that these minuscule modifications 
can add up to measurable time delays for photons from 
cosmological sources. At a redshift z, simultaneously- 
emitted photons, with energies Ei and E2, will arrive at 
the observer with a time delay At = ti — t2 per energy 
difference AE = -Bi - B2 of [9] : 
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where ftm = 0.3, Ha — 0.7 and Hq = 70 km s~^ 
Mpc~^ are the cosmological parameters as currently mea- 
sured. In the case of a vanishing linear term, the mean 
time delay of the photons per squared energy difference 
AB^ =El-El is: 
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The absence of such an energy dispersion has been used 
to set bounds on the parameters ^ and C- Gamma-ray 
bursts and very high energies flares of active galaxies have 
been the primary targets of these "time-of-flight" studies. 
For the linear dispersion term in Eq. ^_^thcsc measure- 
ments reach limits of |C| < 70-150 [13, El, E, El, (3, El 
for gamma ray bursts. For active galaxies, dispersion 
measurements exist for only two sources: Mkn 421 and 
Mkn 501. Both are located at a similar redshift of '^0.03. 
For Mkn 421, a limit of |^| < 200 was set by the Whipple 
collaboration during a flare in 1996 [l3|. For Mkn 501, 
an indication of higher energy photons lagging the lower 
energy ones was reported during a flare in 2005 by the 
MAGIC collaboration 11711 . This dispersion was recently 
quantified to |^| ~ 30 [l^. Since the signal is however 
also marginally consistent with zero dispersion, limits of 
1^1 < 60 and |C| < 2.2 x lO^"^ were derived While hm- 
its in ^ from time-of-flight measurements are approaching 
unity and probing Planck-scale energies, limits on C are 
generally still far from this domain. 

Time-of-flight measurements provide the most direct 
and model- independent test of the constancy of the speed 
of light with energy. However, alternative methods set 
more stringent limits relying on additional assumptions: 
Limits of |^| < 10"'' are deduced if the speed of light 
modifications in Eq. [1] are helicity dependent 1^, 2^ , as 
predicted by some of the models 0, Q . Also, constrain- 
ing limits of If I < 10"^'' and ^ > — 10~^ were recently 
reported in [21| under several assumptions, for example 
the sign of the speed of light modification is assumed to 
be negative or helicity dependent and standard kinemat- 
ics are required to be valid in a Lorentz- violating regime. 

A caveat of time-of-flight measurements is that disper- 
sion might be introduced by intrinsic source effects, which 
could cancel out dispersion due to modifications of the 
speed of light. In the case of a non-detection of disper- 
sion this scenario is unlikely, since it requires both effects 
to have the same time scale and opposite sign. How- 
ever, this "conspiracy of natere" [l6j can only be ruled 
out with certainty by observations of sources at multi- 
ple distances, as - in contrast to dispersion from speed 
of light modifications - source intrinsic dispersion should 
not scale with distance. Population studies of this kind 
have been performed for ganima-ray bursts, resulting in 
limits of If I < 1300 |i,[ll,[3H- For active galaxies the 
data-set is currently too sparse to perform these studies. 
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FIG. 1: Black points show the integral flux VHE light curves 
measured on July 28 from PKS 2155-304 by H.E.S.S. between 
200-800 GeV (upper panel) and >800 GeV (lower panel), 
binned in two-minute time intervals. The zero time point 
is set to MJD 53944.02. Gray points show the oversampled 
light curve, for which the two-minute bins are shifted in units 
of five seconds. The inlay in the upper panel illustrates this 
in a zoom, where the horizontal error bar shows the duration 
of the bin in the original light curve. 



In the present study, photon time delays were searched 
for during the VHE flare of the active galaxy PKS 2155- 
304 observed by the High Energy Stereoscopic System 
(H.E.S.S.) on July 28 in 2006. PKS 2155-304 is located 
at a redshift of z = 0.116 f25|, almost four times more 
distant than Mkn 501 and Mkn 421. The hght curve 
shows fast variability (~ 200 s) and covers an energy 
ran ge o f a few TeV with no significant spectral variabil- 
ity [23| . Considering the unprecedented photon statistics 
10000 photons) at these energies, this flare provides a 
perfect testbed. The data presented here were analyzed 
using the standard H.E.S.S. analysis, described in detail 
in [2J] . Time delays between light curves of different en- 
ergies were sought in order to quantify a possible energy 
dispersion. For this, two different methods were applied, 
which are described in the following. 

The first method determines the time lag between two 
light curves with the Modified Cross Correlation Func- 
tion (MCCF) 01 The MCCF is a standard cross corre- 
lation function [26|, applied to oversampled light curves. 
This allows time delays below the duration of the flux 



bins to be resolved 25|. To optimize the energy gap be- 



tween two energy bands, while keeping good event statis- 
tics in both, the correlation analysis was performed on 
the light curves between 200 and 800 GeV and above 800 
GeV (see Fig. [T]). The mean difference of the photon 
energies between the two bands is 1.0 TeV and the mean 
quadratic difference is 2.0 TeV^. The MCCF of these 
light curves is shown in Fig. [2l In order to measure the 
time delay, the central peak of this distribution was fit- 
ted by a Gaussian function plus a first-degree polynomial, 
resulting in a maximum at Tpeak = 20 s. 




High energy lags- 



:600 ■ 
- 500 ^ 
-400 
: 300 
-200 

: 100 



-600 -400 -200 



200 400 600 
X[S] 



-60 -40 -20 20 40 60 80 100120 



FIG. 2: Left: MCCF of the light curves in Fig. [T] The 
black line shows the best fit of a Gaussian plus first degree 
polynomial. The peak of the fitted function is located at 
Tpcak = 20 s. Right: Cross Correlation Peak Distribution 
(CCPD) obtained from 10000 simulated light curves. The 
shaded area shows the range of the CCPD for Tpcak < 0, 
corresponding to 21% of the total area. The dotted line shows 
the position of Tpeak from the left panel. The CCPD is slightly 
asymmetric, with a mean of 25 s and an RMS of 28 s 



The error on the measured time delay is determined 
by propagating the flux errors via simulations. Ten thou- 
sand simulated light curves were generated for each en- 
ergy band, by varying the flux points of the original over- 
sampled light curve within its measurement errors, tak- 
ing into account bin correlations. For each pair of light 
curves, the peak of the the MCCF was determined, re- 
sulting in a Cross Correlation Peak Distribution (CCPD) 
shown in the right panel of Fig. [2l The CCPD has an 
RMS of 28 s andyields the probability density of the er- 
ror of Tpeak 27|, |28[ . For 21% of the simulations the time 
delay is negative, therefore the measured time delay of 
20 s is not signiflcantly different from zero. 

The response of the MCCF to energy dispersion is 
complex. Primarily, dispersion is expected to shift light 
curves in time according to their mean energy. How- 
ever, dispersion also broadens their structures and pho- 
tons might even get shifted out of a burst, decreasing 
the overall correlation. These "second order" effects 
become increasingly important once the time shifts ap- 
proach the time scale of the observed structures in the 
light curve. The response of the MCCF to dispersion 
was therefore determined by injecting artificial disper- 
sion into the H.E.S.S. data and measuring its effect on 
the CCPD. As shown in Fig. [3l the CCPD follows the in- 
jected time shift per energy linearly in the range of inter- 
est here, confirming the expected behaviour. The second 
order effects mentioned only introduce small deviations, 
visible at higher dispersion values. Nevertheless, the 
measured time delays are transformed to dispersion-per- 
energy with the calibration curve shown in Fig. [3l in or- 
der to take these effects into account. Since the measured 
Tpcak was compatible with zero, a 95% confidence upper 



limit on a linear dispersion of 73 s TeV is given. Ap- 
plying the analogous procedure to a quadratic dispersion 
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FIG. 3: Mean of the CCPD as a function of the dispersion 
injected in the H.E.S.S. data. The points have been shifted 
by the mean value of the CCPD of the original data shown in 
Fig. [2] to emphasize the relative time shifts. Each CCPD is 
derived from ten thousand simulated light curves. The error 
bars show the RMS of distributions. The solid line shows the 
calibration curve used to transform time shifts into dispersion. 
For comparison, the dotted line shows the linear response 
function expected from the mean energy of the correlated light 
curves (see text). 



in energy yields a 95% confidence limit of 41 s TeV ^ . 

The accuracy of the MCCF method was verified with 
an independent set of simulations. Eleven thousand new 
photon lists were generated from the real data using a 
parametric bootstrap method. The parametric model 
was obtained from a polynomial spline fit to the light 
curves in time bins of one minute and a fit of the energy 
distribution of the events in the real data. The CCPD 
of these new simulations confirmed the previously mea- 
sured error on the time delay. Artificially introduced 
dispersion was always recovered within the expected ac- 
curacy. It should be noted that the dispersion limit does 
not depend strongly on the choice of preset parameters, 
such as the energy ranges and time binning of the light 
curves and the fit range of the MCCF peak. Varying 
these parameters within a reasonable range has only a 
small effect (<5 s) on the final result. 

To confirm the result obtained with the MCCF anal- 
ysis, the dispersion measurement was repeated with an 
independent method, which is widely applied in time lag 
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Light curves 



studies of GRB light curves [Ij 
were constructed in two energy bands, and a search for 
extrema was done using a Continuous Wavelet Trans- 
form (CWT) [30]. For this the Last Wave package ^ 
was employed, which provides a list of extrema candi- 
dates with their positions. The extrema were associated 
in pairs between light curves and their relative time delay 
was measured. The association was performed with an 
algorithm based on the Lipschitz coefhcient as in 

The two energy bands were chosen to be 210 to 250 
GeV and above 600 GeV, with a mean energy difference 
of 0.92 TeV. Since tiny dispersions are to be probed, a 



time bin-width of 60 seconds was found to be optimal 
for this study. The CWT method identified two pairs of 
extrema with a mean time delay of 27 seconds. In order 
to assess the error of this value, samples composed of 
hundreds of Monte Carlo experiments were analyzed for 
three linear dispersion values: and ±45 s TeV"^, in 
analogy to the MCCF calibration. The values of the error 
on the measured time lag were found to range between 30 
and 36 seconds. The relation between injected dispersion 
and measured time shift between light curves is again 
used to derive a limit on the dispersion, resulting in a 
95% confidence limit of 100 s TeV^^. The impact of 
systematic effects have also been investigated: selection 
of gamma-like events and the choice of the energy domain 
or time binning of the light curves change the results 
by 0.5tT at most. Various cuts on the CWT parameters 
have been applied and lead to negligible changes in the 
extrema identification. 

The measured limits on the energy dispersion translate 
into limits on the energy scale of speed of light modifi- 
cations. For a linear dispersion in energy, Eq. [2] yields 
1^1 < 17 (or Ep > 7.2 xlO" GeV) for the limit ob- 
tained with the MCCF method, at 95% confidence. The 
linear dispersion limits obtained from the Wavelet anal- 
ysis yields a limit of |^| < 23 (or \^\-^ Ep > 5.2 xlO^"^ 
GeV), confirming this result. These limits are the most 
constraining limits from time-of-flight measurements to 
date. For a quadratic dispersion in energy, the MCCF 
method yields |C| < 7.3 xlO^^ (or ICI"^/^ Ep > 1.4 xlO^ 
GeV) with Eq. EI 

This measurement opens a new redshift range for pop- 
ulation studies of time delays from active galaxies, which 
are needed to rule out the possibility of time delay cancel- 
lation. For a final verdict on this question further VHE 
observations of active galaxies are needed. However, the 
result already shows that the time delay reported for Mkn 
501 in [isll, if considered significant, cannot be attributed 
to speed of light modifications. Current and future in- 
struments such as Fermi for gamma ray bursts, or the 
proposed Cherenkov Telescope Array for active galaxies, 
will further improve the sensitivity of time-of-flight mea- 
surements, perhaps one day revealing deviations from 
Einstein's postulate. 
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